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Abstract. This paper investigates finite-dimensional representations of VT- 
symmetric Hamiltonians. In doing so, it clarifies some of the claims made in earlier 
papers on 'PT-symmetric quantum mechanics. In particular, it is shown here that 
there are two ways to extend real symmetric Hamiltonians into the complex domain: 
0^ , (i) The usual approach is to generalize such Hamiltonians to include complex Hermitian 

Hamiltonians. (ii) Alternatively, one can generalize real symmetric Hamiltonians to 
include complex PT-symmetric Hamiltonians. In the first approach the spectrum 
^ \ remains real, while in the second approach the spectrum remains real if the PT 

symmetry is not broken. Both generalizations give a consistent theory of quantum 
mechanics, but if D > 2, a _D-dimensional Hermitian matrix Hamiltonian has more 
£f) ' arbitrary parameters than a £>-dimensional 'PT-symmetric matrix Hamiltonian. 
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It has been observed that non-Hermitian Hamiltonians that exhibit VT symmetry 
can have real spectra. For example, the class of non-Hermitian Hamiltonians 

H = p 2 + x 2 (tx) u (1) 



CN 



have positive real discrete spectra so long as v > and appropriate boundary conditions 



are specified PjJ El • The domain v > is the region of unbroken VT symmetry, while 
v < is the region of broken VT symmetry. The distinction between these two regions 
is as follows: When v > 0, the eigenstates of H are also eigenstates of VT, but when 
v < 0, the eigenstates of H are not eigenstates of VT . In the unbroken region the 
eigenvalues of H are real and in the broken region some eigenvalues of H may be real, 
but the rest appear as complex-conjugate pairs. 

In a recent letter it was shown that in the region of unbroken VT symmetry a PT- 
symmetric Hamiltonian possesses an additional symmetry represented by the complex 
linear operator C [3]. The operator C commutes with H and with VT and can be used 
to construct an inner product whose associated norm is positive. The theory defined 
by the complex Hamiltonian (JI} with v > is a fully consistent and unitary theory of 
quantum mechanics [3]. 

One might conjecture that VT symmetry is a generalization of Hermiticity. 
However, as we will argue in this paper, this view is not quite precise. Rather, we will 
argue that the appropriate way to construct complex Hamiltonians is to begin with a real 
symmetric Hamiltonian and to extend the matrix elements into the complex domain in 
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such a way that certain constraints are satisfied. There are two distinct ways to perform 
this construction. First, one can generalize real symmetric Hamiltonians to the case of 
Hermitian Hamiltonians and second, one can generalize real symmetric Hamiltonians 
to the case of PT-symmetric Hamiltonians that are not Hermitian. In the second 
generalization the symmetry of the Hamiltonian is maintained but the matrix elements 
are allowed to become complex with the condition that the VT operator commutes with 



Many of the Hermitian Hamiltonians commonly studied in quantum mechanics 
are actually real and symmetric. For example, this is the case of the Hamiltonian 
representing a particle in a real potential V(x), so that H = p 2 + V(x); this Hamiltonian 
is explicitly real 5J. To show that it is symmetric we display it as a continuous matrix 
in coordinate space: 



This matrix is explicitly symmetric under the interchange of x and y. The PT- 
symmetric Hamiltonian in (Q) is also symmetric in coordinate space; however, it is 
complex for all v > 04 

In this paper we investigate the case of finite-dimensional matrix Hamiltonians. We 
show that Hermitian matrix Hamiltonians and PT-symmetric matrix Hamiltonians are 
both acceptable generalizations of real symmetric matrix Hamiltonians. Furthermore, 
they define consistent theories of quantum mechanics. We also demonstrate that for 
the case of D-dimensional matrices the class of Hermitian matrix Hamiltonians is 
much larger than the class of PT-symmetric matrix Hamiltonians. Specifically, we 
know that for large D the number of real parameters in a real symmetric matrix is 
asymptotically \D 2 and the number of real parameters in a Hermitian matrix is D 2 . 
We will see that the number of real parameters in a PT-symmetric matrix Hamiltonian is 
asymptotically \D 2 . The overlap between the classes of Hermitian and PT-symmetric 
matrix Hamiltonians is only the class of real symmetric matrices. A Venn diagram 
showing the relationships between the classes of Hermitian, PT-symmetric, and real 
symmetric matrix Hamiltonians is given in Fig. 

To construct a finite-dimensional PT-symmetric matrix Hamiltonian we begin by 
defining the operators that represent time reversal T and parity P. Both of these 
operators represent discrete reflection symmetries and thus we must have T 2 = P 2 = 1. 
Furthermore, we assume T and P are independent operators, so that they commute 
[P, T] = 0. For simplicity we define the time reversal operator as complex conjugation. 
One can also define T to be Hermitian conjugation (complex conjugation and transpose). 
However, we will see that because all of the relevant matrices in the theory are symmetric 
it makes no difference whether T performs a transpose. It is also possible to choose a 
more complicated definition for T. For example, T could be the combined action of 
complex conjugation and multiplication by some complex matrix. Such alternative 

\ This Hamiltonian is complex even when v is a positive even integer because the boundary conditions 
associated with the eigenvalue problem H<j) = E<f> are complex. See Ref. pQ. 



H. 
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definitions for T will be considered later and we will argue that without any loss of 
generality we may define T to be just complex conjugation. 

Next, we consider the operator V representing parity. Since V commutes with T, 
the entries in the matrix representing V are all real. Furthermore, we will see that V 
must be symmetric. (If it is not symmetric, then we will find that the C matrix that we 
will ultimately construct will not commute with the Hamiltonian H. As a result, the 
quantum theory will violate unitarity. We will return to this point later on.) 

The fact that V 2 = 1 implies that all the eigenvalues of V are either +1 or —1. To 
construct the most general D-dimensional matrix V let us suppose that there are m + 
positive eigenvalues and m_ negative eigenvalues, where m + + m_ = D. That is, if V 
is diagonalized, then it has the form 

P = diag{l,l,..-,1,-1,-1,...,-1}. (3) 

The most general parity matrix can be expressed as 

V = RP Q R~\ (4) 

where R is the most general D- dimensional rotation (orthogonal) matrix [H].§ 

There are \D(D — 1) arbitrary parameters in the matrix R. However, there are 
fewer than this number of parameters in the matrix V in Indeed, if m_ = so 
that Vq is the identity matrix, then there are no arbitrary parameters in V. The exact 

§ Of course, one could take the matrix R to be more general than orthogonal by choosing it to be 
unitary. However, in this case the parity operator V will be complex and will not commute with T. 




Figure 1. Venn diagram showing that the intersection between the classes of 
Hermitian and 'PT-symmetric matrix Hamiltonians is the class of real symmetric 
matrix Hamiltonians. 
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number of arbitrary parameters in V is given by the formula 
-D(D — 1) — -m + {m + — 1) — -m_(m_ — 1). 



(5) 



2 " '2 

Clearly, when D is even, V has the greatest number of arbitrary parameters if 
m + = m_ = |-D. When D is odd, the number of parameters is maximized if we 
choose m + 



m_ 



1; that is, m + = ^(D + 1) and m_ = |(D — 1). Thus, for all D, the 



greatest number of parameters in V is given by the formula 
V-i [l -(-!)»]. 



(6) 



Let us illustrate these results. The most general one- dimensional parity matrix 
V — 1 has no free parameters. The most general two-dimensional parity matrix has one 
parameter: 



V 




(7) 



V 



/ 



(8) 



The most general three-dimensional parity matrix has two parameters: 

' cos 2 — sin 2 cos 29 sin 20 cos 9 — sin 2 sin 29 
V = sin 20 cos 9 — cos 20 sin 20 sin 9 

— sin 2 sin 20 sin 20 sin 9 cos 2 + sin 2 cos 20 

Now, let us construct the most general .D- dimensional PT-symmetric matrix 
Hamiltonian H, where by PT-symmetric we mean that the operator VT commutes 
with H. We will assume that the matrix H is symmetric. (If H were not symmetric, 
then the eigenvectors of H would not be orthogonal. || We will consider the possibility 
of an asymmetric H later.) To count the number of parameters in H we take the parity 
matrix to be in diagonal form Vo as in (jSJ). If the operator VqT commutes with Hq, 



V H* = H n V { 



then Hn has the 2x2 block form 




(9) 



(10) 



where A is a real symmetric m + x m + matrix, C is a real symmetric m_ x m_ 
matrix, and B is a real m + x m_ matrix. Thus, the number of parameters in Hq is 
\D(D + 1). We then transform P to V using the rotation matrix R, and find that the 
number of arbitrary real parameters in the corresponding PT-symmetric Hamiltonian 
H = RH R~ l is given by the combined number \D(D + 1) of free parameters in Hq 
and the number (jHJ) of free parameters in V. Since Ho is not Hermitian (it is complex 
and symmetric) and R is orthogonal, as we have argued above, it follows that the 
Hamiltonian H is non-Hermitian and is not unitarily equivalent to any Hermitian matrix. 

|| The inner product here is just the ordinary dot product, v • v = v T v. 
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Table 1. Number of arbitrary real parameters in the following most general D x D 
matrices: (i) real symmetric parity V, (ii) Po^-symmetric Hq, (iii) 'PT-symmetric H, 
(iv) Hermitian H, and (v) real symmetric H . 



Dimension D 


1 


2 


3 


4 


5 


6 


Large D 


Real Symmetric P: \D 2 - \ [l - (-1)°] 





1 


2 


4 


6 


9 


1 n2 


P T-Symmetric H : \D{D + 1) 


1 


3 


6 


10 


15 


21 


2 


PT-symmetric H: \D 2 + \D - \ 1 - (-1)°] 


1 


4 


8 


14 


21 


30 


3 n2 


Hermitian H: D 2 


1 


4 


9 


16 


25 


36 


D 2 


Real Symmetric H: \D{D + 1) 


1 


3 


6 


10 


15 


21 


2 



As an example, for the case D = 2 the most general PT-symmetric Hamiltonian, 
where P is given in ((7|), contains four free parameters and has the form 

I r + 1 cos — is sin iscos0 + tsin0 \ 
I iscos0 + tsin0 r — tcos0 + is sin / 

The most general 3x3 PT-symmetric Hamiltonian has eight free parameters. 
For arbitrary D there are 

-D 2 + -D--11- (-1)°] (12) 
4 2 g L v ) j v ; 

real parameters in the most general PT-symmetric matrix Hamiltonian. For purposes 
of comparison, in Table Q we give formulas for the number of free parameters in the 
most general real symmetric D x D parity matrix, the most general matrix H Q that 
commutes with VqT, the most general PT-symmetric matrix H, the most general 
Hermitian matrix Hamiltonian, and finally the most general real symmetric matrix 
Hamiltonian. 

Once we have found the most general PT-symmetric matrix Hamiltonian we 
proceed according to the recipe described in Ref. 0]. First, we find the energy 
eigenvalues. The eigenvalues for H in (jll)) are 

e± = r ± tcosa, (13) 

where sin a = s/t and the unbroken PT-symmetric region is s 2 <t 2 . 
Next, we find the corresponding eigenstates: 



. . 1 / sin a cos | — z(l =F cosa) sin ^ \ , ^ , 

a/2(1 =f cos a) cos a I sin a sin \ + z(l =f cos a) cos ^ 

Because we are in the unbroken PT-symmetric region these states are also eigenstates 
of the VT operator. We have chosen the phase in ((T3j) so that the eigenvalue under the 
PT operator is unity: 

VT\s + ) = \s + ), 

TT\s-) = |e_). (15) 
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It seems appropriate now to define an inner product with respect to the VT 
operator. To do so we define the VT conjugate (-| of the state |-) as follows: 

(•I = [PT\.)] T , (16) 

where T is matrix transpose. The VT inner product of two states \a) and |6) is now 
defined as the dot product of the VT conjugate of \a) and \b): 

(a\b) = [VT\a)f.\b). (17) 

This inner product has the symmetry property (a\b)* = (b\a). 

By virtue of (fT5j). for the eigenstates of the Hamiltonian the state (e±\ is just the 
transpose of \e±). The states in (JT5j) are normalized so that their VT norms are 

(e_|e_)=-l. (18) 

Also, the matrix Hamiltonian is symmetric, so these states are orthogonal with respect 
to the VT inner product: 

(e + |e_) = (e_|e + ) = 0. (19) 

Finally, we construct the C operator as outlined in Ref. j3j: 

C= |e+)(e + | + |e_)(e_| 

1 / cos — i sin a sin sin + % sin a cos 
cos a I sin + i sin a cos — cos <p + i sin a sin q. 

It is easy to verify that the matrix C commutes with VT and with H and that C 2 = 1. 
The eigenstates of the Hamiltonian are simultaneously eigenstates of C: 

C\e+) = + \e+), 

C\e.) = - |e_). (21) 

Using these results we can define a new inner product in which the bra states are 
the CVT conjugates of the ket states: 

(•I = [CVT\.)] T . (22) 

The CVT inner product of two states \a) and \b) is now defined as the dot product of 
the CVT conjugate of \a) and |6): 

(a\b) = [CVT\a}} T .\b). (23) 

This inner product has the symmetry property (a\b)* = (b\a). The advantage of the 
CVT inner product is that the associated norm is positive definite. 
We recover the parity operator 



(20) 




that was used in Ref. jl] by choosing = n/2. All the results that are reported in 
Ref. 0| are also obtained for this choice of 0. However, there is an error in Ref. In 
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this reference the parameters s and t in the Hamiltonian must be identical; they cannot 
be unequal because then the matrix would not be symmetric and the eigenvectors would 
not be orthogonal. 

What happens if we choose the parity operator V to have an irregular distribution 
of positive and negative eigenvalues? For example, suppose we take D = 8 and choose 
m + = 6 and m_ = 2. [Of course, in this case there are only twelve real parameters 
in V instead of the sixteen parameters that occur in the symmetric case for which 
m + = m_ = 4. See @. Correspondingly, there are also four fewer parameters in 
the Hamiltonian.] We have found that the signs of the VT norms [see (fT7j) ] of the 
eigenstates of the Hamiltonian are exactly the same as the signs of the eigenvalues of 
V. However, the order of the signs depends on the values of the parameters in H and 
is unpredictable. The operator C is exactly what is needed to cancel each of the minus 
signs in the VT norm so that the CVT norms of the eigenstates are all positive. 

The natural question that arises is whether it is possible to have a more general 
formalism for PT-symmetric matrix Hamiltonians; that is, to have matrix Hamiltonians 
with more arbitrary parameters than the number given in (|T2"j) . There are two 
possibilities: First, one could consider having an asymmetric matrix Hamiltonian H or 
an asymmetric parity matrix V. Second, we could generalize the time reversal operator 
to include a matrix multiplying the complex conjugation operator. 

If the matrix Hamiltonian H is not symmetric, then eigenstates of H corresponding 
to different eigenvalues will not be orthogonal. This forces us to generalize the VT inner 
product (•[•) to include a weight matrix W [Zj. That is, rather than having an ordinary 
dot product of vectors, we would have to generalize the definition of the inner product 
to (-|W|-), where the matrix elements of W are chosen so that 

(e m \W\e n ) = 5 mn . (25) 

In this case, the matrix W plays the same role as the operator C. The drawback of this 
generalization is that W will not commute with H. As we now argue, we must reject this 
generalization of the Hamiltonian because the theory is not unitary: Unitarity means 
that the inner product of two states is independent of time. In the Schrodinger picture 
the states [a, 0) and \b, 0) at time t — evolve into the states \a, t) = e~ lHt \a, 0) and 
\b,t) = e- iHt \b,0) at time t. Thus, 

(a,t\ = [VT\a,t)} T = (a } 0\e iHt . 

The inner product between these states will not be independent of time unless 
e iHtyy e -iHt _ -yy ( rememDer that H commutes with VT), and this requires that W 
and H commute. If W and H do not commute, the theory must be abandoned because 
it violates unitarity and is therefore physically unacceptable. 

Similarly, if we generalize the parity operator to the case of an asymmetric matrix 
V, the most general PT-symmetric H will be asymmetric. Again, we must reject this 
possibility. 

Finally, we ask if it is possible to generalize T so that it is a product of some 
matrix B and the complex conjugation operator. The condition that T 2 = 1 implies 
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that BB* = 1. Also, the requirement that [V, T] = imposes the constraint [P, -B] = 0. 
These two conditions are so strong that no additional parameters appear in the most 
general PT-symmetric matrix Hamiltonian H. 

We do not believe, as has been claimed (see, for example, Ref. jH] and references 
therein), that Hermiticity is a special case of VT symmetry. The problem with the 
analysis in Ref. || is that the norm associated with the inner product is not positive. 
To observe this nonpositivity we construct a vector that is a linear combination of 
eigenvectors of the Hamiltonian: fi\e m ) + v\e n ), where /i and v are complex numbers. 
According to Eq. (15) in Ref. jH], the norm of this vector is /x 2 + z/ 2 , which is not positive 
in general. 

In conclusion, the matrix constructions presented in this paper have changed our 
views regarding the relationship between Hermiticity and VT symmetry. We have 
found that PT-symmetric Hamiltonians should not be regarded as generalizations of 
Hermitian Hamiltonians; rather, based on our study of finite matrices we understand 
that these are two totally distinct and unitarily inequivalent complex classes of 
Hamiltonians whose overlap is restricted to the class of real symmetric Hamiltonians. We 
conjecture that the picture in Fig. Q continues to be valid even for infinite-dimensional 
coordinate-space Hamiltonians. 
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